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Abstract: 



■ Predictions for (</>^)4 theory from renormalization- group- improved perturba- 

^ • tion theory, as formulated by Liischer and Weisz, are compared to published (and 

some unpublished) data from lattice Monte-Carlo simulations of the 4-dimensional 
Ising model. Good agreement is found in all but one respect: — the change in the 
wavefunction-renormalization constant Zji across the phase transition is significantly 
greater than predicted. A related observation is that propagator data in the broken 
phase show deviations from free-propagator form — deviations that become larger, 
not smaller, closer to the continuum limit. More data closer to the critical point are 
needed to clarify the situation. 



1 Introduction 



A complete understanding of the ((/>'*)4 theory is important not only as a funda- 
mental problem in quantum field theory, but also for its implications for the Higgs 
mechanism. According to conventional wisdom the continuum limit of lattice ((/'^)4 
theory is described by a Renormalization Group (RG) analysis using KG functions 
calculated in perturbation theory. The theory was developed in detail by Brezin et 
al ^ and by Liischer and Weisz (LW) |21 01 • It predicts "triviality" in the sense that 
the renormalized coupling gn tends to zero. 

In the late eighties LW's numerical predictions were compared with Monte- 
Carlo data for the 4-dimensional Ising model in both the symmetric phase (^[3] and 
the broken phase [HI • The aim of this paper is to revisit this comparison in the light 
of the much larger data set now available |H1 El El • 

One motivation for this exercise is the recent controversy between Balog, Dun- 
can, Willey, Niedermayer and Weisz (BDWNW) '8_ and Cea, Consoli, and Cosmai 
(CCC) [71. It is important to note that the raw data of the two groups agree very 
well; the dispute is solely over interpretation. I cannot pretend that my sympathies 
are neutral; for many years I have collaborated closely with CCC in gathering lat- 
tice Monte-Carlo evidence (HI ^3 fo'^ unconventional view of "triviality" in (0^)4 
theory advocated by Consoli and myself ^21- Nevertheless, I intend here to take a 
detached view and, except for a few remarks in Sect. 5, I shall not discuss the ideas 
of Ref. 112]. No change in my position is implied; I simply want to focus here on a 
limited question: how well does perturbative RG theory agree with all the available 
lattice data? 

While BDWNW and CCC both describe fits of their data to formulas based 
on 1- or 2-loop perturbation theory, neither group has made a comprehensive com- 
parison to the full LW theory incorporating 3-loop and higher- twist effects. That 
exercise is performed here in Sect. 3, considering all measured quantities in the 
both the broken and symmetric phases. Good agreement is found, except for the 
wavefunction-renormalization constant Zr, as discussed in Sect. 4. Broken-phase 
propagator data are examined in Sect. 5. BDWNWs data show the same deviation 
from free-propagator behaviour found in Ref. 9J. Because of this deviation, differ- 
ent strategies for extracting "mass" and "Z" parameters can lead to very different 
conclusions. Non-Ising data are briefly considered in Sect. 6 and conclusions are 
summarized in Sect. 7. An appendix gives a detailed discussion of the data used. 
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2 Basic definitions 



I shall use the notation of LW and BDWNW and I refer the reader to those papers 
for full definitions. Only a few key facts will be outlined here. The lattice action for 
the (f)'^ theory is written as 



2k ^ ,p{x)<P{x + A) + (Pixf + A {(P{xf - 1) ' 
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LW also define another parameter A that varies between and 1 as A ranges from 
to oo. The limit A — oo (A = 1) corresponds to the Ising model, where <p{x) can 
take only the values ±1. For a given A there is a critical k separating the symmetric 
and broken phases. As k ^ Kc the correlation length (inverse of the physical mass 
in lattice units) diverges, according to the RG theory, so that k ^ Kc corresponds 
to the continuum limit. 

The s/Tk factor between the Ising field <p and the canonical field 00 is the source 
of several notational nuisances. In particular, the field renormalization constant Zr 
defined by LW includes both the trivial 2k factor and the dynamical effects. I 
follow BDWNW in defining = 2kZr as the canonical field renormalization that 
obeys Zr < 1. The renormalized coupling constant ga and renormalized mass 
mR are defined as in LW, as are the susceptibility x and vacuum expectation value 
V = {4>). To plot the data for the latter quantities I have first removed the power-law 
dependence on [k — Kc) by defining 



V^iK 



(4) 



Notice that the combination xv^ is "dimensionless" in this sense. 

The RG theory involves coupled differential equations containing three RG 
functions /3, 7 and S. LW provide a recipe for constructing these functions to 3-loop 
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order, including power-suppressed (or "higher twist") scahng violations to 1-loop 
order. The higher- twist terms, while negligible in the limit important 
further away from Kc- Because of them the integration of the differential equations 
must be done numerically. Using Mathematical I have implemented LW's procedure 
exactly as described in Refs. j2], |SI (referred to below as LW(I) and LW(II), respec- 
tively). ^ The LW procedure depends upon three integration constants Ci,C2,C3 
as well as on the assumed value for Kc- I have verified that my program precisely 
reproduces the results in the LW tables when the same input parameters are used. 
(The Ci constants are defined in the symmetric phase; in the broken phase LW 
define corresponding constants C{, C2, C3 and then prove that C[ = e^'^Ci, while 
C2 = C2 and C3 = C3. I shall quote only numerical values for the Cj's, but of course 
the conversion from Ci to C( is taken into account in my program.) 
In the Ising case a quite precise value for Kc is known |13j : 

Kc = 0.074848(2) (5) 

and will be adopted here. This value is consistent with earlier estimates, 0.074834(15) 
[H] and 0.074851(8) [IS]. [In fact, I initially made fits using 0.074834 and then tried 
0.074851 and found a small but distinct improvement. On closer inspection, the im- 
provement seemed to have slightly "overshot" and I was experimenting with slightly 
smaller "compromise" values when I became aware of the result of Ref. jl3j . Thus, 
I adopt that value both because it has the smallest quoted uncertainty and because 
it seems to produce the best fits.] 

3 Fits to lattice data 

Table 1 of LW(II) gives predictions for the three integration constants based on 
LW(I)'s analysis that matched the RG procedure to the "high-temperature" (small 
K.) expansion in the symmetric phase. For the Ising model these predictions are 

lnCi = 1.5(2), InCa = 1.87(1), InCs = -3.0(1). (6) 

One can indeed fit the available symmetric-phase data quite well with parameters in 
this range (see later). However, these parameter values do not yield a good fit to the 
broken-phase data. This observation is in accord with the experience of CCC [7j. 
^Copies of my Mathematica programs are available upon request. 
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However, BDWNW (HI point out that, with hindsight, the uncertainties quoted by 
LW may have been over-optimistically small, especially in the Ising case. Also, LW 
used a much cruder approximate value for Kc (0.07475(7)) which is another source 
of the discrepancies found by CCC. By adjusting the values of and the Cj's, 
BDWNW claim one can fit the broken-phase data very well with the conventional 
RG theory. To check this assertion, I have considered all available Monte-Carlo 
Ising data (see Appendix) and made a number of fits, adjusting the parameters by 
trial and error. Fig. 1 shows the best fit I have obtained, which uses the following 
parameter values: 

lnCi = 1.24, lnC2 = 1.83, InCg = -2.90. (7) 

Indeed, this fit seems to be an entirely satisfactory description of the broken-phase 
data, given the theoretical uncertainties, especially at the larger k's. 

Comparing these parameter values with the LW values quoted above one sees 
a lower InCi value - in accord with a remark in BDWNW that a value ~ 1.2 is 
needed. However, note also the lower lnC2 value, which will be crucial in what 
follows. 

[It is hard to quote meaningful uncertainties on the "best fit" parameter values 
given above. The effects of varying the Cj's are highly correlated, difficult to describe, 
and often hard to understand intuitively. Also, one should give more weight to fitting 
the data points closer to Kc, where the theoretical uncertainties are less. As a very 
rough guide I would say that changing any of the In Q values in Eq. by plus or 
minus 0.06, 0.01, 0.01, respectively, would lead to a discernible deterioration in the 
quality of the fit and changes by twice these amounts would be unacceptable.] 

However, while the Cj values in Eq. Q yield an excellent fit to the broken- 
phase data, they do not yield a good fit to the symmetric-phase data; see Fig. 2. 
The fits to niR and gn are quite good, but the theoretical curve for Zr lies far below 
the data points. That fact is a direct consequence of the lower C2 value, since Zr 
has a factor of C2. 

If one increases lnC2 to 1.862 (back in accord with LW's predicted value) one 
can then obtain a good fit to symmetric phase data; see Fig. 3. However, this change 
spoils the fit to the broken-phase data, not only for Z^, but also for x and v; see 
Fig. 4. I have tried adjusting the Cj's further but I cannot find any "compromise" 
values that would make the problem go away. One can fit either the broken-phase 
data, or the symmetric-phase data — but one cannot fit both well simultaneously. 
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Figure 1: Broken-phase data for the Ising model compared with the LW procedure 
for parameters InCi = 1.24, lnC2 = 1.83, InCs = -2.90, and Kc = 0.074848. 
The solid curve is the full LW procedure (3-loops plus higher-twist corrections to 
1-loop). The long-dashed lines show the effect of omitting 3- loop terms, while the 
short-dashed lines show the effect of omitting higher-twist terms. The boxed regions 
in (e) and (f) are shown on an expanded scale in (g) and (h). Data points from Refs. 
101 IHl El Uni, see Appendix for details. 
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Figure 2: Symmetric-phase data for the Ising model compared with the LW proce- 
dure for the same parameters as Fig. 1. Note that in (a) the curve hes weh below 
the Zn data points. The boxed region in (c) is shown in more detail in (d). Data 
points from Refs. 01121 El El; see Appendix. 
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Figure 4: Fits to the broken-phase data with the same parameters as in Fig. 3: 
InCi = 1.24, lnC2 = 1.862, InCs = —2.90. Note that these parameters are basicahy 
in accord with LW's predicted values, Eq. ©. Compared to Fig. 1 the fits to 
V, Xv'^, as well as Zj^ in (a)-(d) are spoiled, though there is little or no effect on 
(e)-(h). 
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4 The "step" in i 



The essence of the problem is ihustrated in Figure 5, which shows on both sides 
of the phase transition, combining Figs 2(a) and 1(d) on a common scale. It is 
convenient to define the "step" as: 

A = Zij(K=0.074) — Zr{k=0.075i). (8) 

From the CCCS data point at 0.074 and combining the five data points around 
0.0751 one finds an "experimental" value of A = 0.071(6). However, the theoretical 
curve predicts a step of only about 0.04. Moreover, as argued below, this is a 
robust prediction, essentially independent of the particular Cj values, and with little 
theoretical uncertainty. 




0.0725 0.073 0.0735 0.074 0.0745 0.075 0.0755 0.076 

K 



Figure 5: The "step" in Zji across the phase transition. The data indicate a step of 
0.07, whereas the predicted step is only 0.04. The solid curve is for the parameters 
of Eq. used in Figs. 1 and 2. Increasing C2 and hence Zcrit = 2kcC2 would shift 
the theoretical curve upward with almost no change in shape. See also Fig. 12. 



The theoretical prediction for Zr can be understood simply if we neglect 
higher-twist effects. In the symmetric phase it takes the form 

Zr = 4rit (1^) + 1^" + 0.100896a2 + . . , (9) 

where 

Zcrit = 2kcC2 and a = gij/(167r^). (10) 



9 



At K = 0.074, the measured gn is about 27, so a = 0.171 and the series is well- 
behaved. Numerically one finds, adding on a higher-twist contribution based on the 
fit in Fig. 2: 

Zr{k=o.074) = Zcrit(0.9887)(1.0124) + 0.0080. (11) 

[Note that in Fig. 2(a) the higher- twist contribution is quite sizeable. Curiously, 
it seems to compensate for the linear rise caused by the k/kc factor, so that Zr is 
almost flat until very close to Kc-] In the broken phase one has 

^^ = ^-*(£)(^-^«-°-^^«^^^"' + •••)• ^''^ 
At K = 0.0751 the measured is about 21, so a = 0.133 and again the series is 
well behaved. Numerically, adding on a higher-twist contribution, one obtains 

Zij(«;=o.075i) = Zcrit(1.0034)(0.9646) - 0.0017. (13) 

Hence, the numerical prediction for A is 

A = 0.0331Zcrit + 0.0097. (14) 

The broken-phase data suggest Zcrit ~ 0.933 while the symmetric-phase data imply 
0.963; in any case, Z^rit is certainly less than 1. The higher-twist contribution 
depends in principle on the Cj parameters, but seems to vary little in the various 
fits I have made. There may be some further correction from higher-loop higher-twist 
contributions, but these should be only a fraction of the higher-twist contribution 
already allowed for. Thus, the theoretical prediction for A cannot really be pushed 
above 0.05, well short of the "experimental" value 0.071(6). 

This is not an entirely new problem. It was noted by Jansen et al ,6^ that 
there was a "small discrepancy" for Z/j; their data point at k = 0.076 was about 
2.5(T below the LW prediction. (That is indeed what we see in Fig. 4(d) above.) 
The same problem showed up in Ref. At that time it was not unreasonable to 
suppose that higher-order /higher- twist effects could explain away the discrepancy. 
However, now that there is data much closer to Kc that explanation is no longer very 
credible. 

5 Propagator data 

Another indication that something unconventional may be going on comes from 
data for the momentum-space propagator. Ref. 0, referred to as CCCS below, 
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found that the propagator in the broken phase shows significant deviations from 
free-propagator form. Moreover, those deviations become more evident closer to Kc, 
contrary to conventional ideas about "triviality." In this section I re-examine that 
data and also point out that BDWNW's data show the same effect. 

Some technical preliminaries are needed. is defined as the Fourier trans- 

form of the connected two-point function. I shall normalize to the canonical field 0o; 
hence my G differs by a 2k factor from BDWNW, but agrees with CCCS. On an 
lattice the allowed momenta are = where is a vector with integer- valued 

components and the variable p^, the lattice analogue of the invariant p^p^, is defined 
as 4sin^(p^/2), summed over /x = 0, . . . , 4. 

BDWNW do not directly provide propagator data in the Ising case (though 
they do in two non-Ising cases; see the next section). However, their Table 2 gives 
data for the time-slice correlation function S{t) at k = 0.0751 on a 48^ lattice. Prom 
this data one can construct the momentum-space propagator G{p'^) for momenta 
p = (27r/48)ra (n = 0, 1, 2, . . .) along the time axis by taking the Fourier transform: 

= S{t)e'P* = S{0) + 2 ^ S{t) cos pt + (-1)'^5(24). 

(The property S{t) = 5(48 — t) has been used. The 2k. factor is to convert to the 
canonical normalization.) Without access to the raw data, I am not able to compute 
realistic error bars, but the smoothness of the data, and the good agreement with 
CCCS's data, both suggest that the error bars would be comparable to those of 
CCCS. 

The resulting G{p^) points are plotted in Figure 6 in various ways. BDWNW's 
remark, at the end of Sect. 5, that "the inverse propagator is remarkably linear in 
k"^ up to the maximal (on-axis) momentum = 4"' is seemingly justified by Fig. 
6(b). However, the parameters for this linear fit are quite different from those that 
BDWNW obtained from the low-momentum points (see Fig. 6(a)). 

Deviations from free-propagator behaviour are more easily seen by plotting 
the quantity 

C = C(J5^m) = G(p2)(j52 + m2). (15) 

Of course, (" depends crucially on what one chooses to use as the "mass," m. Fig. 
6(c) uses the = 0.1691 value given by BDWNW in Table 5, while Fig. 6(d) uses 
the larger mass, 0.200, of the fit line in Fig. 6(b) and shows a ^ that is remarkably 
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0.95 



m=0.200 
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Figure 6: Propagator data for k = 0.0751 from BDWNW plotted in various ways. 
No error bars are shown, (a) shows the inverse propagator G^^ with the straight 
hne being BDWNW's free-propagator fit to the first three p ^ points (see Table 
5 of Ref. 0). (b) shows another straight-line fit to the same data that is good at 
almost all momenta. However, the line actually fails to fit the lowest momentum 
points, (c) and (d) show the corresponding ( plots, C = G{p^){p^ + m^), in which 
the deviations from free-propagator behaviour are more easily seen. 
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constant except for a dramatic spike at low momentum — just as seen in Figs. 3, 
4, 5 of CCCS 0. 

It is necessary to distinguish three different "masses." From the conventional 
viewpoint, the physical mass is determined from the exponential fall-off in time of 
the time-slice correlator for zero 3-momentum, S{t). I denote this mass by mTs{0), 
as in CCCS. The measured values of rnTsi^) from CCCS are shown in the second 
column of Table 1. 

In the LW procedure the "renormalized mass" niR and the wavefunction renor- 
malization constant Z/j are defined in terms of an expansion of the inverse propagator 
about p"^ = 0: 

G(p2)-i = Z^i(m| + p2 + 0(p4)). (16) 

This means that C,{p^ ,mR), plotted against p^, should have zero slope at the origin. 
I have used this fact to determine empirically for the CCCS data sets at k = 
0.076, 0.07512, 0.07504, adjusting m until the lowest points lined up. (See Figs. 7, 
8,9.) The resulting ttir values, given in the third column of Table 1, are only slightly 
larger than the mTs{^) values. The ratio between the two is quite consistent with 
the theoretically predicted formula |3J . For most of the subsequent discussion the 
small difference between mTs{^) and mn can be ignored. 

The third "mass," denoted by m-iatt, corresponds, as in Figs. 6(b,d), to the 
mass that gives the best fit to a free-propagator form: 

G{f)-^^Z^,l^{p'' + ml,,), (17) 

at all momenta, excepting the first few low-momentum points. The values found by 
CCCS are given in the last column of Table 1. 



K 


rriTsiO) 


ruR 


mia.u 


0.076 


0.3912(12) 


0.393 


0.42865(456) 


0.07512 


0.1737(24) 


0.176 


0.20623(409) 


0.07504 


0.1419(17) 


0.1426 


0.17229(336) 



Table 1: Measured masses in the broken phase [5]. The zero-momentum time- 
slice mass mTs(O) is the physical mass, from the conventional viewpoint. The rriR 
values, obtained as described in the text, differ slightly by the expected perturbative 
correction. The miatt values are those obtained by CCCS as giving the best fit to 
free-propagator form at all except the very lowest momenta. 
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It is important to note that in the symmetric phase the propagator shows no 
visible deviation from free-field behaviour; see Fig. 1 of CCCS j) for k = 0.074. All 
three versions of the "mass" are indistinguishable. 

However, the situation is quite different in the broken phase. In Figs. 7, 8, 9 
I re-plot CCCS's propagator data using mpt as the mass in Q. With this mass C, has 
zero slope at = 0, by construction. However, C, then rises with j5^; quickly at first, 
then more slowly. The deviation from constancy is highly significant, statistically. 
Moreover, the deviation from free-propagator behaviour is even larger at k = 0.07512 
and 0.07504, closer to Kc, than it is at k = 0.076. 
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Figure 7: Propagator data for k = 0.076 from CCCS plotted as C = G{p^){'p^ + m?') 
for a mass m = mji = 0.393. For this mass ^ has zero slope at = and its value 
there is Zr. 

While the data sets in Figs. 7, 8, 9 are exactly the same as those shown in 
Figs. 3, 4, 5 of CCCS |3, the plots have a completely different appearance. This 
fact is entirely due to the different "m" used in forming just as one sees in Figs. 
6(c), (d) above. In CCCS's figures, which use miatt, the C data points are almost 
exactly constant, except for the lowest 3 or 4 points, which rise up to a dramatic 
peak at = 0. The peak value, C(0) ""^latt), is CCCS's quantity ".^0." 

How this situation comes about is illustrated in Fig. 10, which uses a simple fit 
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Figure 8: As Figure 7 but for k = 0.07512, with mji = 0.176. The scale is exactly 
the same as Figs. 7 and 9. Notice that tlic deviation from constancy is even greater 
than in Fig. 7, even though we are now closer to the continuum limit. 
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0.94 




0.92 



Figure 9: As Figure 7 but for k = 0.07504, with mji = 0.1426. The curve represents 
a simple, empirical fit to the data, to be used in Fig. 10. 
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to the K, = 0.07504 data.^ The lowest curve corresponds to using m = rriTsi^) ~ ''^R^ 
while the uppermost curve corresponds to using m = miattj chosen so that the curve 
is almost exactly flat over the whole range of above, say, 0.1. 




0.1 0.2 0.3 0.4 1 2 3 4 5 6 7 8 




Figure 10: Illustration of the effect of plotting the propagator data as ^ = G{p^){p'^+ 
m^) using different masses m. For clarity, the range of is separated into to 0.4 
and 0.4 to 8. The curves use the empirical fit to the n = 0.07504 data shown in 
the previous figure. The lower solid curve represents C with m = rriTsiO) = 0.1419, 
which differs only slightly from the curve for m = niR = 0.1426 shown in the previous 
figure. The upper solid curve shows C with m = miatt = 0.17229 used in CCCS, Fig. 
5. For this mass C, is essentially constant, as for a free propagator, except for the 
lowest three data points below = 0.1. The dashed curves represent intermediate 
choices of mass m in equal steps of m between mTs(O) = 0.1419 and miatt = 0.17229. 



CCCS argue that their data are indicative of a continuum limit in which, with 
'TT'iatt viewed as the physical mass, the propagator tends to free- field form at all 
finite p^ except for a "spike" in C, at infinitesimally small p^. The top of this spike, 
Zff, = C(0) "T-latt); should diverge to infinity logarithmically in this scenario. Indeed, 
grows from 1.05 to 1.31 between k = 0.076 and 0.07504. The point I want to 
^The fit function corresponds to 

witli parameters A w 0.876, B « 0.081 and M/rriTsiO) close to 3. This form also fits the propagator 
data in the other two cases, with quite similar values of A, B and curiously, the best-fit M is again 
close to 3mTs(0). 
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make here is that the change in viewpoint as to what is the "physical mass" is crucial. 
Ironically, the "odd" features of the data — which, from the conventional viewpoint 
are a too-low Zr associated with a distinct dip in the ^ plots at low momentum — 
are, from the CCCS viewpoint, evidence for a logarithmically growing spike. 

6 Non-Ising data 

BDWNW also collected some data for the theory not in the Ising limit, but 
at A = 0.3 and 0.6, where A is LW's parameter that becomes unity in the Ising 
limit. The other parameters were chosen so that would have a value about 
20, for comparison with their Ising data at n = 0.0751. Their propagator data at 
A = 0.3, 0.6 show no statistically significant deviation from constancy, even if re- 
plotted using (However, their data only extend to = 0.3 and possibly some 
effect might show up if one had data at higher momenta.) 



gR=20 

0.2 0.4 0.6 0.8 1 

X 

Figure 11: Zji in theory at gji = 20, as a function of the parameter A. The 
region between the two dotted curves represents the prediction of LW(II) . The three 
bold points are the data points from BDWNW, Table 5. A deviation is seen only in 
the Ising case, A = 1. 

Figure 11 plots the BDWNW results for in comparison with the LW ex- 
pectation, indicated by the region between the two dotted curves. These curves were 
obtained by finding .Zcrit = 2kcC2 from the C2's of LW(II), Table 1, and the Kc's of 
LW(I), Table 1, and then applying the perturbative correction, 0.967, relating Zr 
to Zcrit at gji = 20. The moral of this plot is that the Zr problem discussed earlier 
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appears to show up only in the Ising case; i.e., it becomes visible only for A above 
0.6. 

7 Summary and Conclusions 

In many respects the RG predictions of LW are impressively successful; they explain 
a large amount of data over a sizeable range. However, on close examination, there 
does appear to be a significant problem: The parameters that fit the broken-phase 
data well (Fig. 1) do not fit the symmetric-phase data for Z/j (Fig. 2(a)). Alterna- 
tively, parameters that fit the symmetric-phase data well (Fig. 3) , and which accord 
well with LW's predicted values, do not fit the broken-phase data (Fig. 4). The core 
of the problem is that the data require a downward step in of 0.07 across the 
phase transition, whereas the theory predicts a step of only about 0.04. This is a 
serious concern because it is the proudest boast of the RG method that it can relate 
the behaviours on each side of the phase transition. 

1 

0.98 
0.96 
. 0.94 
0.92 
0.9 
0.88 
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Figure 12: The "step" in Zr across the phase transition, revisited. The solid curve 
is the RG prediction using the parameters of Figs. 3 and 4 that fit the symmetric- 
phase data well. The open circles represent "Zprop" determined from the C values of 
the large-p^ propagator data. I have used some "editorial license" to omit some data 
points that I consider dubious or uninformative (see discussion in the Appendix). 



Figure 12 represents my best summary of the situation, expressed in conven- 
tional terms. The RG curve here corresponds to Ci parameters chosen to give a 
good fit to the symmetric-phsise data. Relative to this RG prediction, the measured 
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Zrs in the broken phase are too low. The open circles correspond to "^prop" de- 
termined from the large-p^ behaviour of the propagator data. These points appear 
to join smoothly to the symmetric-phase points (where Zr and .^prop a^re indistin- 
guishable). In the broken phase, the vertical gap between the open circles and the 
Zn data points is a measure of the deviation from free-propagator behaviour. Note 
that this gap increases as one approaches Kc- 

Because of the effect illustrated in Fig. 10 these features of the data can be re- 
interpreted, from CCCS's viewpoint, in terms of a logarithmically growing Z^. Prom 
the conventional viewpoint, one can only say that there is a puzzle that remains to 
be resolved. My conclusion is that there is strong motivation for collecting more 
Ising-model lattice data, especially closer to Kc, on both sides, where the residual 
uncertainties from higher-twist effects will be even smaller. 
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Appendix: Data 



The broken-phase data in Fig. 1 come from various sources. Triangles (a) represent 
data for x and ^ from CCC fT. Stars ("^) represent data from BDWNW |Sj and also 
an earlier data point from Jansen a/ jHj at k = 0.076. A convenient compilation of 
these data points can be found in Table 3 of BDWNW [HI. Diamonds (♦) represent 
data at K = 0.07512 and 0.07504 for Z^, g^^ and m/j that I have extracted from 
results of CCCS Ref. 0. (The CCCS results at 0.076 agree completely with Jansen 
et al.) 

The symmetric-phase data in Fig. 2 also come from various sources. Triangles 
(a) represent data from Montvay, Miinster, and Wolff Stars (■^) represent earlier 
data from Montvay and Weisz The diamond (♦) at k = 0.074 comes from CCCS 
The diamond at k = 0.0724 represents unpublished data of Cea and Cosmai 
|10j : see comments below. 

In general there appears to be very satisfactory agreement between the data 
of the various groups. (To avoid clutter I have generally not plotted a data point 
that completely agrees with, but is less precise than, the equivalent data point from 
another group.) A few data points, however, deserve comment because superficially 
they might appear to weaken my case for a large step in Zr. 

The BDWNW data point for at k = 0.0754 appears high in comparison 
with the others. I believe that the most likely explanation is an unlucky 2a statistical 
fluctuation. I say this (and it implies no criticism of BDWNW) for two reasons: (i) 
Zn in the region 0.075 < n < 0.076 is expected to vary slowly and smoothly, so at 
least one of the BDWNW points must be shifted by more than l.ba. Since the other 
BDWNW points are well corroborated by independent data, the 0.0754 point is the 
likely culprit, (ii) The 0.0754 data point also seems to lie slightly off the fit curves 
in the m/j and qr plots ((e) and (f) in either Fig. 1 or Fig. 4); the hypothesis that 
the rriR value has a 2a upward fluctuation would explain away all the discrepancies. 

In fact, after having written the previous paragraph, I learned of new data 
of Cea and Cosmai ^Oj at the same k, which indeed yield a smaller mass, = 
0.262(1) and a lower Zr = 0.898(8). This result is included in the Zr figures as a 
square (■). (To avoid clutter it has not been included in the mn plots.) 

In the symmetric phase, the k = 0.0724 result of Montvay et al [Sj for Zr 
appears to be considerably lower than its two neighbouring points, also from Mont- 
vay et al. The quoted errors are very small, so statistics should not be a factor. 
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The predicted Zn in this region is almost exactly constant, so the sharp dip implied 
by the three Montvay et al (a) points (see Fig 2(a) or 3(a)), if it were real, would 
be in serious disagreement with theory. Because of this anomaly I asked Cea and 
Cosmai to repeat the Monte-Carlo calculation of Z/j at the three k values studied by 
Montvay et al. They kindly did so |Tn] and found excellent agreement at 0.0710 and 
0.0732, but at 0.0724 they found a different result that indeed interpolates smoothly 
between the neighbouring k's. I therefore suspect that there must be some trivial 
mistake in Montvay et aFs 0.0724 point - perhaps in the transcription of the actual 
computer result to the published paper. 

The symmetric-phase data point at k = 0.074, closest to the phase transition, 
is obviously important. The Montvay- Weisz (MW) data point (■^) for mji comes 
from the last line of Table 2a of W as 0.2125(10) times a factor 1.002 (see last 
sentence of Sect. 3.1) to convert the physical (or "time-slice") mass to m/j. This 
gives 0.2129(10), which agrees very weh with CCCS's result 0.2141(28) (not shown 
in the figures). The MW result for Z/j also comes from the last line of Table 2a of 
MW, from "z = 6.44(10)" converted by a 2k factor and (1.002)^ (see Eqs. (26), (27) 
in MW), giving Zji = 0.957(15). This is compatible with the more precise value 
from CCCS, 0.9682(23) obtained from propagator data. The gn value from MW 
comes from remarks in their Sect. 4.2 that a gji of 24 with a 10 — 15% possible 
variation would fit their data. 
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